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摘 要 : 为 了 更 好 的 表现 自然 图 像 的 曲线 特性 ， 进 一 步 提 高 数字 水 印 算法 的 鲁 棒 性 ， 提 出 一 种 基于 Curvelet-DSVD 和 
视觉 密码 相 结 合 的 强 鲁 棒 零 水 印 算法 。 首 先 对 原始 图 像 进 行 Arnold 置 乱 ; 其 次 进行 Curvelet 变换 得 到 低频 域 信息 ， 对 
低频 域 信息 进行 分 块 并 对 各 个 块 进行 双 奇 异 值 分 解 (DSVD), 利用 块 最 大 奇异 值 与 整体 奇异 值 均值 之 间 的 关系 构造 特征 
答 阵 ， 同 时 利用 视觉 密码 将 水 印信 息 生 成 两 个 共享 份 ; 最 后 将 其 中 一 个 共享 份 进行 Arnold 置 乱 后 与 特征 矩阵 进行 异 或 
运算 生成 替 水 印 。 实 验 结果 表明 ， 该 算法 能 够 有 效 地 抵抗 常规 攻击 ， 与 现 有 的 零 水 印 算 法 相 比 ， 和 鲁 棒 性 更 强 ， 安 全 性 
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Strong robust zero watermarking algorithm based on Curvelet - DSVD and visual cryptography 


Qu Changbo, Wu Deyang 
(College of Software Liaoning Technical University, HuLudao Liaoning 125105, China) 


Abstract: In order to better perform the characteristics of natural image curves and further improve the robustness of digital 
watermarking algorithms, this paper proposed a robust zero-watermarking algorithm based on Curvelet-DSVD and visual 
cryptography. Firstly, performing Arnold scrambling on the original image, then obtaining the low frequency domain information 
by Curvelet transform, segmenting the low frequency domain information, and each block was double singular value 
decomposition(DSVD) , by using the relationship between the maximum singular value of the block and the mean value of the 
global singular value structuring characteristic matrix. At the same time, this paper used the visual cryptography to generate two 
shared copies of the watermark information. Finally, one of the shares was execute to Arnold scramble and then xor operation 
with the feature matrix to generate a zero-watermark. Experimental results show that the proposed algorithm can effectively 
resist the conventional attacks and is more robust and more secure than the existing zero-watermarking algorithms. 
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抵抗 几何 攻击 的 影响 ， 但 滤波 攻击 和 噪声 攻击 对 构造 的 特征 和 矩 
阵 产 生 一 定 影响 。 文献 [6] 提 出 一 种 基于 DWT-SVD 的 图 像 双 零 

为 了 解决 传统 变换 域 水 印 算 法 的 透明 性 与 鲁 棒 性 之 间 的 矛 。 水 印 算法 ， 利 用 haar 小 波 对 载体 图 像 进行 两 次 变换 ， 选 择 载体 
盾 ， 温 泉 等 人 提出 零 水 印 的 概念 ， 由 于 零 水 印 算法 不 需要 将 。 图像 的 低频 域 LL2 来 构造 零 水 印 , 随 着 对 水 印 算法 的 不 断 研究 ， 
版 权 信息 嵌入 到 载体 图 像 中 ， 对 载体 图 像 的 内 容 完整 性 起 到 很 。 近年 来 出 现 许 多 较为 新 颖 的 水 印 算法 TI， 文 献 [7] 选 取 图 像 感 
好 的 保护 作用 ， 所 以 近年 来 零 水 印 在 版 权 保护 领域 得 到 广泛 的 ”兴趣 的 区 域 构造 特征 矩阵 ， 再 将 特征 抢 阵 与 水 印 图 像 进行 异 或 
应 用 己 3]。 运算 生成 零 水 印 。 算 法 对 于 常规 攻击 具有 一 定 的 鲁 棒 性 ， 但 当 
传统 的 零 水 印 算法 大 都 利用 小 波 变 换 和 奇异 值 分 解构 造 载 。” 感 兴趣 域 受 到 攻击 时 ,水 印信 息影 响 较 大 ,具有 一 定 的 局 限 性 。 
体 图 像 的 特征 信息 , 然后 与 水 印信 息 结合 生成 零 水 印 H9, 文献 ”文献 [8] 利 用 余弦 离散 小 波 变换 (DCT)〉 具 有 能 量 集中 优势 来 构 
[各 利用 离散 小 波 和 奇异 值 分 解构 造 零 水 印 , 对 于 小 范围 的 噪声 。” 造 零 水 印 ， 但 利用 图 像 块 均值 构造 特征 矩阵 稳定 性 较 差 。 为 提 
攻击 、 滤 波 攻击 和 压缩 攻击 具有 较 好 的 稳定 性 ， 但 对 于 较 大 强 ”高 图 像 特征 的 稳定 性 ， 曾 文 权 等 人 外 利用 整数 小 波 变 换 和 混沌 
度 的 攻击 鲁 棒 性 较 差 。 文 献 [5] 提 出 一 种 基于 不 变质 心 的 鲁 棒 零 加密 原理 提出 一 种 基于 整数 小 波 变 换 的 鲁 棒 零 水 印 算法 。 实 验 
水 印 算法 ， 选 择 不 变质 心 作为 稳定 的 几何 参考 点 ， 能 够 有 效 的 ”结果 表明 算法 对 于 非 几何 攻击 具有 较 强 的 鲁 棒 性 ， 但 由 于 算法 
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未 考虑 几何 攻击 对 图 像 的 影响 ， 所 以 在 受到 剪 切 攻 击 和 旋转 攻 
击 时 ， 提 取 的 水 印 相似 度 较 低 。 为 了 克服 几何 攻击 对 载体 图 像 
的 影响 ， 曲 长 波 等 人 0 利用 小 波 变换 结合 视觉 密码 提出 一 种 基 
于 视觉 密码 和 边缘 检测 的 零 水 算法 。 算 法 对 于 小 范围 的 噪声 攻 
击 具 有 较 好 的 鲁 棒 性 ， 但 对 于 几何 攻击 表现 鲁 棒 性 较 差 。 肖 振 
和 久 等 人 0 提出 了 一 种 基于 增强 奇异 值 分 解 和 细胞 神经 网 络 的 
零 水 印 算法 ， 算 法 通过 均匀 化 奇异 值 并 结合 神经 网 络 算法 解决 
了 对 角 线 失真 问题 ， 但 神经 网 络 的 引入 ， 增 加 了 算法 的 时 间 复 
杂 度 。 以 上 算法 虽然 使 用 小 波 变换 能 够 很 好 地 将 图 像 分 为 高 频 
或 和 低频 域 ， 但 是 自然 图 像 中 包含 有 大 量 的 纹理 特征 ,曲线 奇 
异性 表现 比较 突出 ,小 波 变换 对 于 点 奇异 特性 表现 较 好 , 而 对 于 
像 曲线 奇异 特性 表现 较 差 ， 所 以 在 受到 攻击 时 提取 的 水 印信 
息 往往 影响 比较 大 。 

针对 上 述 问题 ， 本 文 提 出 一 种 基于 Curvelet-DSVD 和 视觉 
密码 的 强 鲁 棒 零 水 印 算法 。 首 先 对 载体 图 像 进行 Amold 置 乱 ; 
其 次 对 置 乱 后 的 载体 图 像 进行 Curvelet 变换 提取 出 图 像 曲线 特 
征 信 息 ， 并 对 得 到 的 曲线 信息 的 进行 分 块 和 DSVD; 然后 利用 


互 
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1 是 Curvelet 变换 的 分 解 实例 ， 其 中 图 1(a) 大 小 为 512 
X512 的 Lena 图 像 ，(b) 为 Curvelet 变换 分 解 的 图 像 ，(c) 为 
小 波 变换 分 解 的 图 像 。 


6 个 尺度 


图 1 Curvelet 变换 和 小 波 变 换 实例 


图 1 可 以 看 到 , Lena 图 像 经 过 Curvelet 变换 后 被 划分 为 


层 ,第 一 层 为 Coarse 尺度 层 ,第 六 层 即 最 外 层 为 Fine 尺 


页 系数 组 成 的 矩阵 。 


高 由 


度 层 ， 中 间 的 第 二 至 第 五 层 为 Detail 尺度 层 ， 每 层 的 系数 分 别 
被 划分 为 8、8、16、16 个 小 方向 , 都 是 由 


每 层 的 系数 的 绝对 值 的 平方 和 作为 能 和 


而 Lena 图 像 经 过 小 波 变 换 则 被 分 成 4 个 层 ， 有 具体 如 图 1(c) 所 


示 。 分 别 统 计 Curvelet 变换 和 小 波 变换 后 每 层 系数 的 能 量 ， 将 


进行 直方 图 对 比 ， 有 具体 


块 最 大 奇异 值 与 块 整体 最 大 奇异 值 的 均值 之 间 的 关系 构造 特征 
和 矩阵， 同时 对 水 印 图 像 进行 视觉 密码 加 密 操 作 ， 产 生 两 个 秘密 
图 份 , 将 其 中 一 个 秘密 图 份 进行 Amold 置 乱 ， 将 特征 矩阵 与 置 
乱 后 的 秘密 图 份 进行 异 或 操作 生成 零 水 印 ， 最 后 将 生成 的 零 水 
印 与 另 一 秘密 图 份 保存 至 版 权 保护 中 心 。 


1 ”基础 理论 


1.1 ”Curvelet 变换 理论 
Curvelet 变换 是 一 种 多 分 辨 、 带 通 、 具 有 方向 性 的 函数 分 
析 方 法 ， 是 通过 一 种 特殊 的 滤波 过 程 和 多 尺度 Ridgelet 变换 来 


如 图 2 所 示 。 
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(a)Curvelet 变换 能 ] 


其 中 图 
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方 图 (b) 小 波 变 换 能 量 直方 图 


| 
3 4 


2 


图 2 Curvelet 变换 和 小 波 变换 能 量 


直方 图 对 比 


2(a)(b) 横 向 表示 频 域 变换 后 的 层 数 ， 


三 | 


实现 的 ， 对 于 图 像 曲线 边缘 的 描述 ， 其 逼近 性 能 具有 近乎 最 优 
的 非 线 性 允 近 误差 衰减 阶 ， 对 曲线 特征 信息 具有 很 好 的 表示 性 
能 0 。 以 笛 卡 尔 坐标 系 下 的 fl4,41O0<4,6 < 四 作为 输入 ,Curvelet 
变换 的 离散 形式 为 


c2 (lk) = > Flt,t Ibilt,t] (1) 
Ozh rn <n 


其 中 : je042,3.…n 为 尺度 参数 ;1e0.12,3…n 为 方向 参数 ; 
k=(k,k)eZ? 是 一 个 旋转 参数 。Canddes 和 Donoho 提出 了 一 种 
基于 USFFT(unequally spacefast Fourier trans form) 的 快速 离散 
Curvelet 变换 实现 方法 03， 步 又 如 下 : 

a) 对 于 给 定 一 个 笛 卡尔 坐标 下 的 二 维 函数 进行 2DFFT 变 
换 ， 得 到 二 维 频 域 信息 。 


fln, mm /127% 12 (2) 


b) 对 得 到 的 二 维 频 域 信息 进行 重 采 样 , 得 到 相应 的 采样 值 。 


flmsm nm,tan Ql, mm) ep G3) 
c) 将 内 插 后 的 7 与 六 函数 相 乘 得 到 fn,n] 
flms1]= fmt —n tan IU, ln, na] (4) 


对 ,进行 下 FT 逆 变 换 得 到 离散 Curvelet 变换 的 系数 集合 


Cc? (7 。 


图 2 的 能 量 分 布 


纵向 表示 对 应 的 能 


可 知 ， 经 Curvelet 变换 后 的 图 像 的 系 


色 
主要 集中 于 第 一 层 ， 其 余 的 各 层 分 布 的 能 量 较 少 ， 而 经 


过 小 波 变 换 后 的 图 像 的 能 量 分 布 2、3 和 4 层 (主要 含有 高 频 信 


息 ) 的 


1.2” 双 奇异 值 分 解 “DSVD) 
SVD 和 矩阵 分 析 中 常用 的 数学 工具 , 它 可 以 将 一 个 矩阵 分 解 
为 三 个 矩阵 的 乘积 03]， 对 于 一 个 大 小 为 m 


值 分 解 可 表示 为 
6 
2 
T=U,、 
61 
和， 
a 
5 
其 中 : UV 左 奇异 矩阵 ， 
0 


和 > 吃 0 =, >0 3 Vi 为 右 奇 异 矩 阵 。 


Xn 的 矩阵 7 ， 奇 异 


为 奇异 值 矩 阵 ; 


让 


2 


图 像 经 过 奇异 值 分 


解 会 把 元 余 的 信息 去 掉 , 用 少量 的 奇异 值 代 表 整 个 图 像 的 信息 ， 


后 的 图 


在 图 像 去 噪 和 图 


像 压 缩 上 共有 重要 的 意义 ， 同 时 经 过 奇异 值 分 解 


像 具 有 一 定 抗 攻击 能 力 ， 所 以 经 常用 


于 数字 水 印 领域 来 
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提高 算法 的 鲁 棒 性 ， 但 是 传统 奇异 值 分 解 的 安全 性 较 低 且 存在 
一 定 虚 警 问题 0。 

针对 传统 奇异 分 解 存在 的 问题 ， 提 出 一 种 具有 安全 性 更 高 


的 双 奇 异 值 分 解 (DSVD) 。DSVD 不 直接 对 图 像 进行 SVD 分 


解 ， 而 且 先 对 图 像 进 行 双 对 角 化 操作 ， 然 后 对 双 对 角 甜 阵 进行 
奇异 值 分 解 。DSVD 分 解 过 程 如 下 : 
a) 先 对 矩阵 ,进行 双 对 角 化 分 解 ， 如 式 (6) 所 示 。 

Tusa = U: ZI Vi (6) 
其 中 : UV 为 正 交 矩阵，WVW 为 西 矩 阵 。 

i 

pe 才 
下 人 (7) 
天 


Din 
b) 然 后 对 得 到 的 双 对 角 和 矩阵 进行 奇异 值 分 解 ， 得 到 最 大 
奇异 值 矩阵 ; 


ms (®) 
1.3 视觉 密码 
视觉 密码 是 两 位 伟大 的 数学 家 Moni Naor 和 Adi Shamirt!5] 
提出 的 一 种 视觉 系统 解密 的 新 型 密码 方案 ， 它 的 核心 思想 是 ， 
将 一 个 秘密 图 像 (secretimage)，, 利用 视觉 密码 加 密 规 则 加 密 ， 
得 到 两 张 共享 图 份 ， 其 中 第 一 张 共 享 份 可 以 看 做 是 秘密 图 像 加 
密 后 得 到 的 密 文 (cipher) ,第 二 张 共享 图 份 可 以 看 做 是 解码 的 
秘 钥 (key) 。 由 于 每 一 个 共享 份 都 是 像 随 机 噪声 一 样 杂乱 无 章 


的 ， 所 以 密码 破译 者 是 无 法 从 一 张 图 份 中 得 到 另 一 张 共享 份 的 
任何 信息 的 ， 更 加 无 法 破解 秘密 图 像 。 


本 文 利用 的 视觉 密码 原理 将 水 印信 息 产 生 两 个 密 钥 ， 将 其 
中 一 个 密 钥 进行 Arnold 置 乱 后 与 载体 图 像 特征 生成 零 水 印 , 另 
个 密 钥 用 于 后 期 验证 。 其 加 密 过 程 如 式 〈9) 一 〈10) 所 示 。 
1 0| [fo 1 
| 下 | ' | keyl 
1 1| 


| M 


| | 人 


[1 0| [0 1| key2 
[Io 1| |1 0 


(10) 


有 


: 两 个 密 钥 


Im 


否则 为 白色 像素 。 图 3 区 实例 ， 其 中 图 3(a) 为 水 
印 图 像 , (b) 和 (ce) 分 别 是 通过 视觉 密码 原理 生成 的 两 个 互 不 重 羞 
的 密 钥 sharel 和 share2 。 


(c)share2 


(a) 版 权 水 印 (b)sharel 


各 得 到 | ， | , 则 原始 图 像 的 像素 为 黑色 像素 


全 | 
曲 长 波 ， 等 : 基于 Curvelet- nsvcbin 8 全 候 时 法 
2 零 水 印 算法 
2.1 零 水 印 生成 
零 水 印 的 生成 过 程 需 要 构造 图 像 的 特征 信息 ， 首 先 对 大 小 


为 512 X512 图 像 进行 Arnold 置 乱 ， 然 后 对 置 乱 后 的 图 像 进 行 
Curvelet 变换 、 分 块 和 DSVD, 最 后 构造 特征 和 矩阵。 同时 对 


版 权 


水 印 进行 视觉 密码 加 密生 成 两 个 共享 份 ， 将 其 中 一 个 共享 


份 进 


行 Arnold 置 乱 , 最 后 将 置 乱 后 的 共享 份 与 特征 矩阵 进行 异 
作 生 成 零 水 印 ， 将 
中 心 。 零 水 印 具 体 生 成 步骤 如 下 所 示 : 

a) 对 大 小 为 512X512 的 载体 图 像 7 进行 Amold 置 乱 
乱 方式 如 式 (11) 所 示 。 


[all 


mod(N) 


bp) 对 
变换 的 曲线 系数 矩阵 clysewzsc 。 
c 然 后 按照 式 (12) 对 C1256x256 进行 分 块 » 得 


IFmod(i,4) =—=0&&mod(j,4)=0 Then 
B,,{i/4,j/4}= 


cl(m,n) 


其 中 : 3):)),ne{(j -3):(D)} ; 

qd) 根据 式 (6) 和 (8) 求 得 各 个 块 B. 的 最 大 奇异 值 5,， 
最 大 奇异 值 矩阵 $Y ,然后 根据 式 (13) 求 得 最 大 奇异 值 和 矩阵 
值 WM 。 


ijel,2,3...256; mel{(i— kel,2,3:...64 。 


或 操 


生成 的 零 水 印 和 另 一 共享 份 保存 至 版 权 保护 


» 


(11) 


乱 后 的 载体 图 像 进行 Curvelet 变换 ， 得 到 Curvelet 


到 块 矩 阵 B 。 


(12) 


构成 
的 均 


(13) 


e) 利 
的 关系 构造 特征 矩阵 7 。 


} 
yO 
. 


IF SP > M 
ELSE 


每 个 块 最 大 奇异 值 %7 与 整体 最 大 奇异 值 的 均值 4 


(14) 


人 利用 视觉 密码 原理 对 水 印 图 像 W 进行 加 密 操 作 ， 产 
个 秘密 图 份 ww 和 w,， 将 其 中 一 秘密 图 份 w 进 行 Arnold 置 
得 到 置 乱 后 的 秘密 图 份 w 。 


将 人 中 置 乱 后 的 秘密 图 份 mw 与 日 中 得 到 的 特征 矩阵 7 进 


行 异 或 操作 生成 零 水 印 zZ， 并 将 生成 的 零 水 印 z 和 秘密 图 
存放 于 版 权 保护 中 心 ， 两 个 置 乱 密 钥 K, 和 过 由 客户 保管 。 
Z = XORCw ,7) 
零 水 印 产 生 过 程 如 图 4 所 示 。 
2.2 版 权 认证 
版 权 认 证 过 程 需 要 从 版 权 保护 中 心 提取 零 水 印 z 和 秘密 


份 ,版 权 认 证 过 程 的 曲线 特征 信息 提取 过 程 与 2.1 节 的 
a) 一 e) 一 致 得 到 新 的 特征 和 矩阵 了 7。 具体 步 又 如 下 : 


a) 从 版 权 保护 中 心得 到 零 水 印 z 与 新 构造 的 特征 矩阵 
了 异 或 操作 得 到 置 乱 后 秘密 图 份 w 。 
wi = XOR(Z,T') 
b) 从 客户 中 得 到 密 钥 &: 将 秘密 图 份 w 进行 Arnold 逆 
得 到 解密 后 的 秘密 图 份 w 。 


生 两 
乱 ， 


步 又 


7 进 


(16) 
乱 ， 
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c) 将 解密 后 的 秘密 图 份 wi 与 版 权 中 心 的 w, 进行 重 又 ， 得 到 
版 权 水 印 w 。 


载体 图 像 版 权 水 印 
| | 
Arnold 置 乱 | 密 钥 K1 视觉 密码 
VY VY 
了 Share1 Share2 
Curvelet 变 换 


9 © 
分 块 处 理 并 进 


行 DSVD 版 权 水 印 | 版 权 保护 中 心 
Y 
生成 特征 矩阵 Arnold 置 乱 < 一 密 钥 K2 


一 一 >| 异 或 操作 | 


| 


零 水 印 
图 4 零 水 印 生成 过 程 


版 权 认证 过 程 如 图 5 所 示 。 


受到 攻击 的 
载体 图 像 


v 
Arnold 置 乱 ， 一 一 密 钥 K1 


v 
Curvelet 变 换 


分 磊 处 还 并 乒 权 保护 中 
进行 DSVD 人心 

] | 
特征 矩阵 零 水 印 

一 ~ 异 或 操作 到 

Arnold 逆 置 乱 “| 一 密 钥 K2 
! 
Share1 Share2 


— 


版 权 水 印 


图 5 版 权 认证 过 程 


3 ”参数 说 明 及 鲁 棒 性 攻击 实验 


3.1 参数 说 明 

为 了 验证 本 文 算法 的 有 效 性 , 选择 在 64 位 Windows7 操作 
系统 和 MATLAB R2014a 平台 上 进行 仿真 实验 ， 实 验 载体 图 
像 选择 大 小 为 512X512 的 Airplane、Lena、Baboon 、Paper 和 
Bride 的 标准 灰 度 图 像 ， 分 别 对 应 图 6 (9 一 〈e) ;版 权 水 印 图 
像 选择 大 小 为 64x64 含有 ”版 权 水 印 ” 四 个 字 的 二 值 水 印 图 像 ， 
对 应 图 6 (从; 同时 利用 视觉 密码 将 版 权 水 印 图 像 加 密生 成 两 个 
t 享 图 份 ， 分别 对 应 图 6 (g) 和 (h); 为 了 进一步 提高 算法 的 安全 
性 和 和 鲁 棒 性 ， 将 得 到 的 sharel 进一步 置 乱 ， 用 于 增强 算法 的 鲁 


4 6(， 


(d) (©) (9) 


(g)sharel (h)share2 Q) 置 乱 后 的 sharel 
图 6 载体 图 像 与 版 权 水 印 
本 文 用 归 一 化 相关 系数 (NC) 衡量 提取 的 水 印 与 原始 水 印 
相似 度 。NC 是 衡量 数字 水 印 的 相关 标准 , 本 文采 用 如 下 计算 方 
式 [16], 
YY wl 
NC(Ww,w) = = 一 一 一 一 (17) 
志方 之 2 让 


其 中 : w、w 分 别 表示 初始 水 印 和 提取 的 水 印 。 
3.2 ” 虚 警 率 实验 

为 了 验证 算法 的 虚 警 问题 ， 计 算 图 6 的 五 幅 载 体 图 像 构造 
的 零 水 印 之 间 的 相似 性 ， 不 同 零 水 印 之 间 相 似 度 越 低 ， 表 明 虚 
警 率 越 低 ， 否 则 相反 。 有 具体 实验 结果 如 表 1 所 示 。 


表 1 虚 警 实验 结果 (NC 值 ) 


偏 移 参 数 ”Airplane Lena Baboon Paper Bride 
Airplane 1.0000 0.5777 0.4952 0.4849 0.5667 
Lena 0.5777 1.0000 0.4852 0.4480 0.4385 
Baboon 0.4952 0.4852 1.0000 0.5337 0.5184 
Paper 0.4849 0.4480 0.5337 1.0000 0.4847 
Bride 0.5667 0.4385 0.5184 0.4847 1.0000 


表 1 可知， 五 幅 不 同 载 体 图 像 之 间 生 成 的 零 水 印 的 相似 
度 都 很 低 ， 最 高 相似 度 为 0.577 7， 最 低 相 似 度 为 0.438 5, 平均 
相似 度 为 0.503 3， 表 明 本 文 算法 虚 警 较 低 。 
3.3 非 几 何 鲁 棒 攻 击 实验 

常见 的 非 几 何 鲁 棒 性 攻击 主要 为 噪声 攻击 、 滤 波 攻击 和 压 
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缩 攻击 等 。 本 文 实验 主要 硬 


攻击 实验 。 


1) 噪声 攻击 


绕 以 上 几 种 非 几何 攻击 进行 鲁 棒 性 


本 文 对 五 幅 不 同 的 载体 图 像 进行 的 不 同 强度 的 椒盐 噪声 攻 
击 检验 算法 对 于 不 同 强度 噪声 的 抵抗 能 力 ， 并 计算 提取 水 印 的 
NC 值 。 具 体 实验 结果 如 表 2 所 示 。 
表 2 不 同 强度 椒盐 噪声 攻击 实验 数据 表 (NC 值 ) 
噪声 强度 Airplane Lena Baboon Paper Bride 
0.1 0.9751 0.9769 0.95087 0.9838 0.9709 
0.2 0.9531 0.9612 0.9463 0.9549 0.9518 
0.3 0.9308 0.9436 0.9102 0.9399 0.9296 
0.4 0.9288 0.9103 ”0.9038 0.9198 0.9086 
0.5 0.8846 0.8938 0.8601 0.89335 0.8865 


表 2 可 以 看 出 ， 随 着 噪声 强度 的 增 大 ， 从 五 幅 载 体 图像 


中 提取 的 水 印 与 原始 版 权 水 印 的 相似 度 (NC 值 ) 都 很 高 , 平均 


值 在 0.910 0 以 上 。 对 于 载体 图 


像 Lena 对 于 以 上 各 种 强度 的 椒 


盐 噪 声 攻击 时 得 到 的 NC 值 最 高 为 0.9769， 最 低 为 0.893 8。 而 
对 于 载体 图 像 Airplane、Baboon、Paper 和 Bride 得 到 的 NC 值 
的 平均 值 也 在 0.912 3 以 上 。 对 于 以 上 几 种 不 同 强度 的 椒盐 噪 


声 攻 击 ， 整 体 表现 出 较 好 的 鲁 棒 性 能 。 

2)JPEG 压缩 攻击 

对 五 幅 载体 图 像 进行 不 同 强度 的 JPEG 压缩 攻击 ， 压 缩 攻 
击 主要 去 除 载 体 图像 的 元 余 的 信息 ， 载 体 图 像 在 受到 压缩 攻击 


后 ， 生 成 的 零 水 印 也 因此 受到 较 大 的 影响 。 具 体 实验 数据 如 表 


3 所 示 。 由 表 3 


喜 


J] 知 ， 对 于 不 同 的 


压缩 强度 攻击 ， 从 五 幅 载 体 
图 像 中 提取 的 水 印 与 原水 印 的 相似 度 的 平均 值 为 0.994 5 以 上 ， 


表明 本 文 算 法 对 不 同 强度 的 了 PEG 压缩 攻击 表现 出 较 强 的 鲁 棒 
性 能 。 
表 3 不 同 强度 JPEG 压缩 攻击 实验 数据 表 (NC 值 ) 
压缩 因 寺 Airplane Lena Baboon Paper Bride 
10 0.9909 0.9922 0.9778 0.9935 0.9904 
20 0.9939 0.9939 0.9839 0.9974 0.9918 
30 0.9961 0.9974 0.9874 0.9974 0.9935 
40 0.9970 0.9961 0.9848 0.9983 0.9926 
50 0.9957 0.9961 0.9918 0.9987 0.9944 
60 0.9974 0.9987 0.9927 0.9987 0.9939 
3) 中 值 滤波 攻击 
对 五 幅 载 体 图 像 分 别 进行 中 值 滤波 攻击 ， 选 择 的 模板 大 小 


分 别 为 3X3、5X5、7X7 和 9X9 的 滤波 模板 进行 


具体 实验 数据 如 表 4 所 示 。 
不 同 的 滤波 攻击 后 ， 得 到 NC 
X9 的 滤波 模板 得 到 的 NC 值 


的 了 


仿真 实验 。 


表 4 可 知 ， 五 由 
值 在 0.963 7 一 0.999 1， 即 使 是 9 
均值 也 在 0.9813 以 上 。 所 以 


载体 


图 像 在 受到 


本 文 算 法 对 以 上 四 种 不 同 模板 大 小 的 中 值 滤波 攻击 具有 较 强 的 


抵抗 攻击 能 力 。 
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表 4 中 值 滤波 攻击 实验 数据 表 (NC 值 ) 
模板 大 小 Airplane Lena Baboon Paper Bride 
3x3 0.9965 0.9974 0.9852 0.9991 0.9926 
Sx5 0.9917 0.9948 0.9769 0.9970 0.9865 
7x7 0.9891 0.9926 0.9708 0.9957 0.9791 
9x9 0.9856 0.9909 0.9637 0.9913 0.9751 


3.4 ”几何 鲁 棒 攻击 实验 


几何 攻 了 


载体 图 像 发 4 


二 主要 为 旋转 攻击 和 行 
很 大 的 变化 ， 所 以 对 提 
1) 旋 转 攻击 


实验 对 五 幅 载 体 图 


度 分 别 为 1"、 
看 出 ， 本 文 算法 对 于 旋转 攻击 表现 出 较 好 的 鲁 棒 必 


2° 和 3°。 具 体 实 


像 分 别 进行 小 


列 


局 度 的 旋转 攻击 ， 旋 转 f 


移 攻 击 ， 几何 攻击 使 得 
取 的 水 印 影响 也 较 大 。 


验 数据 如 表 5 所 示 。 由 表 5 可 以 


表 5 旋转 攻击 实验 数据 表 (NC 值 ) 


台 
月 EC 。 


旋转 度数 ”Airplane 


Lena 


Baboon 


Paper 


Bride 


向 左 1" 
向 右 1" 
向 左 2" 
向 右 2" 
向 左 3 
向 右 3 


0.9850 
0.97981 
0.9417 ”0.951 
0.9465 
0.9300 


0.9321 


0.9838 
0.9738 


0.9423 
0.9369 
0.9281 


0.9816 
0.9838 
6 0.9466 
0.9537 
0.9376 
0.9325 


0.9833 
0.9887 
0.9427 
0.9403 
0.9280 
0.9363 


0.9816 
0.9724 
0.9416 
0.9430 
0.9288 
0.9348 


2) 行列 偏 移 攻击 
对 五 幅 载 体 图 像 分 别 进 


同行 列 的 移动 攻击 ,移动 丢失 的 像素 进行 


行 上 、 下 、 左 、 


据 如 表 


最 低 为 


多 方向 性 ,月 


6 所 示 。 1 


取 的 水 印 NC 值 在 0.9671 以 上 。 而 
0.923 6, 最 高 为 0.9$06。 由 于 Curvelet 变换 
以 对 于 各 个 方向 的 攻击 时 ， 依 然 能 
水 印信 息 ， 说 明 本 文 算法 对 于 抵抗 各 个 方向 的 行 允 


现 出 的 鲁 棒 性 较 强 ， 


表 6 可 知 ， 对 五 幅 载 体 图 像 人 


表 6 行列 侦 移 攻 了 


篇 移 5 行 时 ， 


扁 移 2 
得 到 的 NC 值 


右 四 个 方向 进行 不 
0 操作 。 


体 实验 数 
列 时 ， 提 


-起 


多 尺度 、 


Ff 实验 结果 (NC 值 ) 


很 好 地 提取 


1 偏 移 攻击 表 


偏 移 
方向 


Airplane 


0.9673 
| -0.9326 
| -0.9716 
0.9260 
0.9718 
0.9251 
0.9846 
0.9331 


Lena 


0.9756 
0.9496 
0.9723 
0.9482 
0.9792 
0.9430 
0.9809 
0.9506 


Baboon 


0.9707 
0.9239 
0.9671 
0.9226 
0.9711 
0.9236 
0.9712 
0.9341 


Paper 


0.9698 
0.9243 
0.9673 
0.9206 
0.9835 
0.9496 
0.9738 
0.9459 


Bride 


0.9638 
0.9416 
0.9695 
0.9361 
0.9726 
0.9420 
0.9752 
0.9413 


3.5 


组 合 攻 击 


为 了 进 


步 验 证 算法 的 强 鲁 棒 怕 


E 能 ， 本 文 增加 组 合 攻 击 仿 
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真实 验 ， 组 合 攻 了 
缩 攻击 + 噪声 攻击 、JPEG 压缩 攻击 + 中 值 滤波 攻击 、 


fF 的 类 型 分 别 为 旋转 攻击 + 噪声 攻击 、JPEG 压 


行列 偏 移 


攻击 + 中 值 滤波 攻击 和 行列 仿 


移 攻 击 + 高 斯 噪声 攻击 等 。 


验 数据 如 表 7 所 示 。 
波 ， 载 体 图 像 Airplane、Lena、Paper 和 Bride 得 到 的 NC 值 都 
在 0.991 3 以 上 。 对 于 JPEG 


表 7 可 知 ， 对 于 JPE 压缩 30%+ 吕 


玉 缩 30%b+ 高 斯 噪声 0.05 攻 玫 


体 实 
P 值 滤 


二 时 ， 


整体 得 到 的 NC 值 的 习 .9668。 总 体 而 言 ， 本 文 算法 对 于 
以 上 几 种 组 合 攻击 表 玫 的 鲁 棒 性 能 。 
表 7 组 合 攻击 实验 数据 表 (NC 值 ) 
攻击 类 型 Lena Baboon Paper Bride 
向 左旋 转 2°+ 椒 
0.9312 0.9261 0.9248 0.9213 
盐 噪声 0.05 
JPEG 压缩 30%+ 
的 0.9756 0.9611 0.9709 0.9582 
JPEG 压缩 30%+ 
0.9970 0.9822 0.9956 0.9913 
中 值 滤波 3x3 
句 下 偏 移 2 行 + 
0.9792 0.9606 0.9839 0.9766 
中 值 滤波 5x5 
句 右 偏 移 2 列 + 
0.9556 0.9931 0.9561 0.9526 


高 斯 噪声 0.05 


3.6 ”对 比 实验 


为 了 更 好 地 验 订 


E 本 文 算法 的 鲁 棒 性 能 ， 选 择 大 小 为 512 X 
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C 合作 上 
曲 长 波 ， 等 : 基于 Curvelet-DSVD 和 视觉 寺 训 的 强 专 棒 替 水 印 算 法 


移 攻击 ， 相 比 文献 [4,8,9] 提 升 %2-9%。 文 献 [4] 使 | 


j 块 最 大 奇异 


512 的 Lena 图 像 作为 载体 图 像 ， 水 印 图 像 选择 大 小 为 64 X64 


含有 “版 权 水 印 ” 的 二 值 氏 
不 同 攻 击 类 型 与 实验 结果 妇 
1) 鲁 棒 性 分 析 

表 8 可 知 ， 本 文 算法 对 于 抵抗 非 几 何 攻 击 和 几何 攻击 相 
能 更 优秀 。 对 于 非 几何 攻击 , 当 


比 于 文献 [4,8,9] 表 现 昌 


像 与 文献 [4,8,9] 进 行 实验 对 比 ， 具体 
[0 表 8 所 示 。 


噪声 强度 为 0.2 时 ， 文 献 [ 和 得 到 的 NC 值 为 0.8040， 重 棒 性 较 


差 。 椒盐 噪声 医 


由 图 7 可 以 看 出 ， 随 着 噪声 强 


度 越 来 越 大 , 本 文 算法 的 NC 值 一 直 比 文献 [4,8,9] 的 要 高 .JPEG 


压缩 攻击 如 图 8 所 示 。 从 医 


8 可 以 看 出 ， 在 不 同 的 JPEG 压 


缩 


强度 下 ， 文 献 [8,9] 与 本 文 算法 得 到 的 NC 值 几乎 一 致 ， 而 对 比 


文献 [和 具有 较 大 的 提升 。1 
局 部 局 限 性 ， 不 具有 方向 性 


六 献 [4,8,9] 使 用 的 小 波 变 换 具有 
， 在 受到 非 几何 攻击 时 ， 图 像 边 缘 


信息 常常 被 削弱 。 本 文 算法 所 使 用 的 Curvelet 变换 和 DSVD 有 


效 地 增强 了 算法 的 鲁 棒 性 
起 到 了 增强 的 效果 ， 相 
更 强 的 表现 能 力 ， 所 L 
对 于 儿 何 攻击 , 
切 攻 击 和 旋转 攻击 方 
升 6% 左 右 。 旋转 攻 


为 Curvelet 变换 对 图 像 边缘 信息 
于 小 波 变换 对 于 图 像 的 边缘 信息 具有 


到 像 在 受到 非 几 何 攻击 时 鲁 棒 性 更 强 。 
文献 [4,8,9] 具 有 较 大 的 提升 , 尤其 在 剪 


面 ， 在 剪 切 攻击 方面 本 文 算法 比 文献 [9] 提 


吧 9 所 示 。 从 图 9 可 以 看 出 , 对 于 不 同 
角度 的 旋转 攻击 , 本 文 算法 的 NC 值 一 直 比 文献 [4,8,9] 的 要 高 ， 


平均 高 出 5% 左右 ， 而 向 左旋 转 3* 时 ， 本 文 的 NC 值 在 0.93 左 


右 ， 文献 [4,8,9] 的 NC 值 者 


了 0.9 以下， 和 鲁 棒 性 较 差 。 对 于 行列 偏 


值 构成 图 像 的 特征 ， 在 像素 位 置 发 生变 化 时 ， 块 的 最 大 奇异 值 
也 会 产生 较 大 的 变化 ,文献 [8] 使 用 DCT 的 均值 构造 的 零 水 印 ， 
在 受到 剪 切 操作 时 ， 图 像 块 均值 影响 较 大 。 而 文献 [9] 直 接 在 空 
或 中 利用 块 均 值 与 整体 均值 的 关系 构造 特征 和 矩阵， 在 受到 几何 
攻击 时 ， 图 像 的 像素 或 像素 的 位 置 发 生 了 变化 ， 直 接 影响 了 块 
均值 与 整体 均值 的 关系 ， 稳 定性 较 差 。 本 文 使 用 的 块 最 大 奇异 
值 与 整体 最 大 奇异 值 构造 特征 ， 最 大 奇异 值 具有 能 量 集中 稳定 
性 能 ,， 相 比 文献 [4,.8,9] 具 有 稳定 性 更 高 。 总 体 而 言 , 本文 算法 对 
于 常规 的 鲁 棒 性 攻击 要 比 文献 [4,8,9] 的 更 具有 优势 ， 鲁 棒 性 更 
强 。 
表 8 本 文 算法 与 文献 [4,8,9] 对 比 实验 数据 表 (NC 值 ) 
攻击 文献 文献 本 文 
参数 文献 [8] 
类 型 [4] [9] 算法 
椒盐 0.05 0.9140 0.9738 0.9739 0.9878 
噪声 0.2 0.8062 ”0.9478 0.9365 0.9612 
高 斯 0.05 0.9454 0.9585 0.9949 0.9843 
噪声 0.2 0.8040 0.9144 0.9707 0.9631 
中 值 
模板 大 小 3x3 0.9819 0.9766 0.9971 0.9974 
高 斯 
模板 大 小 3x3 0.9832 0.9782 0.9949 0.9965 
JPEG 玉 缩 因子 10 0.9300 0.9917 0.9895 0.9922 
玉 缩 因子 30 0.9707 0.9913 0.9927 0.9935 
左上 角 1/16 0.9587 0.9357 0.8712 0.9389 
右 下 角 116 0.9590 ”0.9295 0.8812 0.9300 
向 左旋 转 3° 0.8893 0.8772 0.8474 0.9369 
向 右 旋转 3° 0.8818 0.86677 0.8662 0.9281 
下 移 2 行 0.9033 0.9597 0.9668 0.9809 
右 移 2 列 0.8803 0.9516 0.9597 0.9723 
先 缩小 0.5 倍 再 放 
0.9924 0.9984 0.9984 0.9988 
大 2 倍 
先 放大 4 倍 再 缩 
0.9982 0.9991 0.9993 0.9995 
小 0.25 倍 
1- 
本 文 
;文献 [和 
小 > 文献 [8] 
0.95 * 文献 [9] 
9 0.9 
081 02 04 05 06 
噪声 强度 
图 7 椒盐 噪声 攻击 
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图 9 旋转 攻 了 


It 


在 安 


全 改 


[4,8,9] 只 使 用 传统 


征 过 程 中 


和 


生 , 二 是 加 强 算法 


FE 能 方面 ， 本 文 算法 相 比 文 


的 加 密 操 作 , 安全 1 


对 载体 图 像 进 


天 操作 ， 


时 ， 在 不 知道 特征 提取 的 密 铀 | 


言 息 。 同 时 使 用 视觉 


密 处 理 ， 


密码 和 Arnold 置 乱 方式 对 水 E 


大 [4,8,9] 更 安全 ， 文 献 
生 较 低 。 本文 算法 在 提取 特 
一 是 消除 图 像 像 素 间 的 相关 
的 安全 性 ,攻击 者 在 对 载体 图 像 进行 特征 


提取 


青 况 下 是 很 难 提取 到 


生成 的 零 水 印 


不 是 直接 由 版 权 水 印 构造 的 ， 


用 视觉 密 


码 产 生 的 


似 噪声 点 
中 的 版 权 


吏 


的 密 钥 


信息 。 


份 ， 


4 ”结束 语 


本 文 


基于 Curvelet-DSVD 逢 
克服 了 传统 零 水 印 算 法 不 能 表征 


mm 


基于 Curvelet 变换 、 
I 视觉 密码 的 强 鲁 棒 零 水 印 算法 ， 


图 份 进行 
使 获得 零 水 印信 息 也 无 法 


下 


k 享 图 份 是 一 


[以 即 


DSVD 和 视觉 密码 理论 所 


自然 图 像 


强度 攻击 鲁 棒 性 差 的 问题 。 
分 为 二 ， 实 现 
进行 置 乱 ， 进 


三 方 认证 方式 ，Arnold 变换 对 
步 提 高 了 水 印 的 安全 性 。 
法 对 于 常规 的 噪声 攻击 、 滤 ; 


实验 结果 表 昌 


列 偏 移 攻击 表现 出 较 强 的 鲁 棒 性 能 ， 


同时 算法 的 安全 性 


致 的 特征 
信息 进行 加 


而 是 利 
张 类 
得 知 其 


出 一 种 
该 算法 


线 特征 和 对 于 较 大 
同时 ， 视 觉 密码 实现 将 版 权 水 印 一 
视觉 密码 的 共享 份 


NN 


该 算 


攻击 、 压 缩 攻击 、 旋 转 攻 击 、 行 


E 较 高 。 


chinaxiy 会 作 期 刊 


曲 长 波 ， 等 : 基于 Curvelet-DSVD 和 视觉 密码 的 强 便 棒 零 水 


但 本 文 算法 对 于 较 大 强度 的 组 合 攻击 时 表现 的 鲁 棒 性 稍 差 ， 有 


待 进一步 的 优化 ， 在 下 一 步 的 研究 中 ， 将 进一步 提高 算法 对 于 

复杂 环境 中 的 抗 攻击 能 力 。 
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